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The spectrum of primordial fluctuations from inflation can be obtained using a mathematically 
controlled, and systematically extendable, uniform approximation. Closed-form expressions for 
power spectra and spectral indices may be found without making explicit slow-roll assumptions. 
Here we provide details of our previous calculations, extend the results beyond leading order in the 
approximation, and derive general error bounds for power spectra and spectral indices. Already at 
next-to- leading order, the errors in calculating the power spectrum are less than a per cent. This 
meets the accuracy requirement for interpreting next-generation CMB observations. 
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I. INTRODUCTION 

Recent observations of the cosmic microwave back- 
ground (CMB) P, 01 and of the galaxy distribution in 
real and redshift space 'S' have dramatically improved 
our knowledge of the large-scale Universe. Results from 
these observations have been largely consistent with the 
inflationary paradigm of cosmology, since it appears that 
the Universe is at critical density and that the observed 
structure in the Universe arose from the gravitational col- 
lapse of adiabatic, Gaussian, and nearly-scale-invariant 
primordial density fluctuations, a key prediction of sim- 
ple inflation models. 

The idea that an inflationary expansion, aside from 
solving cosmological puzzles such as the horizon and flat- 
ness problems, also provides a natural mechanism for the 
generation of primordial fluctuations was realized around 
the time the model was first proposed [J . Early estimates 
already established the essential adiabatic, Gaussian, and 
nearly-scale invariant nature of the perturbation spec- 
trum. However, these estimates were rather crude and 
their true value lay in providing very significant quali- 
tative guidance rather than accurate quantitative infor- 
mation. As the prospect of testing inflation and even 
distinguishing between individual models from precision 
CMB observations became more realistic, attention was 
increasingly focused on improving the understanding and 
predictive control of the theory underlying the generation 
of fluctuations. The development of the gauge invariant 
treatment of cosmological perturbations [a, B 13) Si and 
its application to infiation ^ ^ Jli, ilj was of great help 
in clarifying and systematizing the calculations involved. 

The simple nature of the perturbation spectrum char- 
acteristic of infiation arose from models where infiation 
was caused by the dynamics of a single scalar field in a 
relatively featureless potential, evolving in an effectively 
friction-dominated 'slow-roll' regime 12, 13J. It was re- 
alized that in more complex inflation models (multiple 
fields, nontrivial potentials) it is possible to engineer the 



spectra in a variety of ways including violation of scale- 
invariance and introduction of isocurvature fiuctuations. 
However, data from present-day observations do not de- 
mand the consideration of dynamically exotic inflation 
models. 

Inflation predicts a spectrum of metric perturbations in 
the scalar (density) and tensor (gravitational wave) sec- 
tors, the vector component being naturally suppressed. 
Both scalar and tensor perturbations cause anisotropy 
in the CMB temperature |lj]. Scalar fiuctuations seed 
formation of large-scale structure in the Universe, while 
the tensor perturbations lead to a stochastic gravitational 
wave background 15]. In addition, tensor modes cause a 
potentially observable polarization of the CMB, a target 
for next-generation CMB observations [l3- Given the 
variety, precision, and volume of data now and soon to 
be available from the various probes of primordial fluc- 
tuations, it is important to provide controlled theoretical 
analyses for various classes of inflation models. 

One obvious application of such an analysis lies in com- 
paring theoretical predictions to observations in order to 
test specific infiation models. In addition, there is the 
prospect of obtaining information about the slow-roll po- 
tential via a controlled inverse analysis of observational 
data, the so-called 'reconstruction' program [l3| (for a 
review see (Ij)- Finally, there is also the more general 
question of whether the infiationary paradigm itself can 
be tested from observations. Regarding this question, re- 
cent activity has focused on identifying classes of inflation 
models distinguished by their differing observational sig- 
natures [l8j rather than trying to reconstruct a speciflc 
inflationary potential. A recent discussion of parameter 
estimation using input from inflation calculations can be 
found in Ref. ^. 

To expand further on the above, we note that while 
inflation is the most compelling cosmological paradigm, 
it is not without its competitors 20] . None of the alter- 
native scenarios, including structure formation via topo- 
logical defects (definitively ruled out by observations) , or 



string-inspired models such as the ekpyrotic and cychc 
Universe scenarios, as yet may be viewed as offering con- 
vincing competition. Nevertheless, an important ques- 
tion arises in the comparison of predictions particular to 
generic models of inflation to predictions one might ex- 
pect on purely general grounds or from other models. 
A good example of this is provided by the spectral in- 
dex ns characterizing the scalar perturbation spectrum. 
The Harrison-Zeldovich choice ng = 1 [2l| is aestheti- 
cally natural, independent of inflation. A high-precision 
observational detection of ng provides a good target for 
inflation models if ng is indeed slightly different from 
unity. (Clearly, an intuitive argument such as that of 
Harrison and Zeldovich cannot be refined further with- 
out a corresponding theoretical framework.) Related to 
this, another important point is that inflation predictions 
for the scalar and tensor spectra are not mutually inde- 
pendent. Future measurements of the CMB can provide 
tests of this dependency, thus it is important to state 
precisely the appropriate 'consistency relations' [23I and 
the accuracy to which they can be calculated for a range 
of inflation models. 

The equations governing the evolution of the scalar 
and tensor perturbations can be solved numerically. For 
specific models, it can scarcely be argued that an approx- 
imate analytical result is competitive with a mode-by- 
mode numerical solution of the equations. Nevertheless 
analytic results can be extremely powerful in providing 
generic results valid for many classes of models. Accurate 
analytic results also greatly reduce the computational 
requirements for the reconstruction program. There- 
fore some effort has concentrated on improving theoret- 
ical predictions by sharpening the slow-roll analysis. In 
this regard, approximations based on certain slow-roll as- 
sumptions have been criticized on multiple fronts. It has 
been argued that error control is not adequate and that 
they are not systematically improvable [23, |2J|. There 
have been some efforts to alleviate this situation J2^ |2g , 
however, they are based on restrictive assumptions (e.g., 
ns ~ 1) and often lead to complicated mathematical for- 
mulations. 

Recently, we presented an analysis of the inflationary 
perturbation spectrum for single-field models based on 
applying a uniform approximation to the relevant mode 
equations [23. This analysis goes beyond the slow-roll 
approximation, does not make restrictive assumptions, 
has definite error bounds, and is systematically improv- 
able. In this paper we implement the method beyond 
leading order and provide general error bounds on power 
spectra and spectral indices. While our general expres- 
sions for the spectral indices are nonlocal, we can derive 
local expressions for ns and ut by employing a further 
approximation. We also demonstrate that by Taylor- 
expansion of our local results, expressions of the slow- 
roll type can easily be obtained. Finally, in order to 
demonstrate the accuracy of the uniform approximation 
explicitly, we discuss its application to an exactly solvable 
inflation model. 



The organization of the paper is as follows: Section UTI 
provides the essential background to understanding the 
problem of computing the power spectrum. Section IIIII 
explains the uniform approximation, Section llVl gives the 
leading order results, Section provides the results at 
next order, and Section IVII contains the results of Sec- 
tion IIVI simplified to a local form. Here we also encap- 
sulate the conventional slow-roll results and derive the 
analogous results from the uniform approximation. Sec- 
tion [^H] considers a special example which possesses an 
exact solution and uses this to demonstrate the charac- 
teristic features of the uniform approximation. In Sec- 
tion ^^H] we comment briefiy on Bessel function approx- 
imants. Section Hxl concludes with a discussion of the re- 
sults obtained and an outline of future work. Appendix IXI 
provides a short list of the relevant definitions and equa- 
tions for linearized perturbations in an FRW universe. 
Appendix IbI contains the technical details behind the er- 
ror formulae of the uniform approximation. 



II. BACKGROUND 

The cosmic microwave background and the large- 
scale distribution of matter both provide information on 
the primordial perturbation spectrum, albeit 'processed' 
by physics during the radiation and matter-dominated 
phases. Extraction of information from large scale struc- 
ture observations is complicated by the nonlinear effects 
of the gravitational instability and various sources of ob- 
servational bias. Fortunately, the physics of the photon 
distribution is very different from that of the mass distri- 
bution: As a consequence of radiation pressure, the CMB 
is very uniform and fluctuations in it can be treated ad- 
equately by a linearized analysis. 

The generation of perturbations during inflation is due 
to the amplification of quantum vacuum fiuctuations by 
the dynamics of the background spacetime. Of course, 
the post-infiationary epoch must also be understood in 
order to make a connection with observations. A serious 
treatment of the dynamics of the post-infiationary phase 
is necessarily quite complicated as various effects such 
as reheating, particle decay, etc. have to be taken into 
account. Fortunately, it turns out that an accurate ac- 
counting of much of this physics is not required in order 
to make inflation predictions for CMB anisotropics and 
the large-scale distribution of matter. 

The modern treatment of fluctuations generated by in- 
flation is in turn based on the gauge-invariant treatment 
of linearized fluctuations in the metric and fleld quan- 
tities (Appendix 13 [SIESH)- A particularly conve- 
nient quantity for characterizing the perturbations is the 
intrinsic curvature perturbation of the comoving hyper- 
surfaces [ll|, C, = ujz^ where u = a{54> -\- cjy'A/h) and 
z = a4>' /h, with a, the scale factor, the prime denoting a 
derivative with respect to conformal time, h = a' /a, 54> 
the perturbation in the homogeneous background scalar 
field, (j), and A^ a quantity characterizing a perturba- 



tion of the background metric (details are given in Ap- 
pendix^) . As shown in Appendix^ the gauge- invariant 
quantity, m, satisfies the dynamical equation 



u — Au u = 0. 

z 



(1) 



It follows immediately that C is approximately constant 
in the long wavelength limit A; ^ 0. This is true during 
the inflationary phase as well as in the post-reheating 
era. Moreover, the Einstein equations can be used to 
connect the gravitational potential $a and ( so that a 
computation of the power spectrum of C provides all the 
information needed (aside from the transfer functions) to 
extract the temperature anisotropy of the CMB. Details 
of this procedure can be found in Refs. [a, |2J, |23. We 
will return to some aspects of this analysis below. 

The calculation of the relevant power spectra involves 
a computation of the two-point functions for the appro- 
priate quantum operators, e.g., 



(0|w(77,x)u(77,x + r)|0) = 
the operator u being written as 



dk sin kr 
k kr 
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where ak, a^ are annihilation and creation operators, re- 
spectively, such that [dk, dj,,] = Skk' , and defining the vac- 
uum state dk\0) = Vfc. The complex amplitude Uk{ri) 
satisfies 



z 

z 



Uk = 0. 



(4) 



Solving Eqn. Q) is the fundamental problem in deter- 
mining the primordial power spectrum Pu (or Pq). The 
corresponding mode equation for tensor perturbations is 
given by 



a 
a 



Vk = Q. 



(5) 



Equations iQJ and jSJ have the mathematical form of 
Schrodinger equations. A simple approach to analyti- 
cal approximation of Eqns. Q and ([SJl relies on the 
fact that exact solutions exist in the limits k^ ^ \z" / z\, 
|a"/a| (short wavelength) and k^ ^ \z" / z\, \a" /a\ (long 
wavelength) or, as will be made more explicit below, as 
—krj — > oo and kr] ^ 0~ . For scalar perturbations, 



Uk 



Uk 



2k 



-*'='' (fc2»|z"/^|, -kr, 



Akz {e<.\z"lz\). 



oo), (6) 
(7) 



Here, the short wavelength solution corresponds to the 
choice of an adiabatic vacuum for modes on length scales 
much smaller than the scale set by the curvature. The 
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FIG. 1: Sketch of the potential barrier for density perturba- 
tions. The vertical dashed lines delineate the three different 
regions where the solution for Uk is investigated as explained 
in the text. 77 marks the point where the fc-mode enters into 
the potential barrier, the turning point. Inflation lasts over a 
wide 77-range, extending Region II as indicated by the arrows. 
In the far right of Region III radiation dominates. 



long wavelength solutions correspond to the growing 
mode on scales much larger than the Hubble length. 
(Analogous solutions exist for gravitational wave pertur- 
bations.) 

The long wavelength solution for C^k — Uk/ z is just the 
(/c-dependent) constant Ak of Eqn. (jT)). In order to deter- 
mine the corresponding power spectrum P(^{k) = Ps{k), 
the main task is to fix the unknown constant Ak by con- 
necting the two asymptotic solutions. We sketch the sit- 
uation in Fig. ^ The large-Zc regime lies in region I and 
the small- A; regime lies in region HI. The two asymp- 
totic solutions may be connected either by a matching 
procedure performed in region II, typical of the slow-roll 
class of approximations discussed further below, or, as 
performed here, by constructing a global interpolating 
solution. Once Ak is determined, the power spectrum 
for C can be computed in the regime {k — > 0) of actual 
interest. 
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where the index S indicates scalar perturbations. 

III. THE UNIFORM APPROXIMATION 

We begin our discussion of the uniform approximation 
by making the substitution 



z" 1 
z rj^ 



in Eqn. Q, yielding 
u'M - 



C\r,) 



Ukiv) = 0. 



(9) 



(10) 



A similar substitution can be made for the case of tensor 
perturbations. 

If, in Eqn. I|1U|I . we make the assumption that C is con- 
stant, (which is the case for power law inflation, see Sec- 
tion VII) then an exact solution in terms of Bessel func- 
tions is immediate. The problem is to solve the equation 
when C is not constant but is a slowly varying function 
of time. Our aim is to do this without being forced to 
state C{r]) in any real detail. The method we will use 
is the technique of uniform approximation as presented 
by Olver jJQ] ■ We now provide a brief summary of this 
treatment for the differential equation of interest H10() . 

The differential equation we wish to solve is of the 
general form 



^^[b^9iv)+qiri)] 



(11) 



Depending on the behavior of h^g{if) +q{'i]) this equation 
has different approximating solutions. In the case that 
g{fj) = at the point fj, so ^ is a turning point, the solu- 
tion is expressed in terms of Airy functions; if 5(77) has a 
pole of order n> 2 the Liouville- Green (LG) approxima- 
tion must be employed ,29(j. The Liouville- Green approx- 
imation is essentially equivalent to the familiar Wentzel- 
Kramers-Brillouin (WKB) approximation; for an appli- 
cation of WKB to computing perturbations from infla- 
tion, see Ref. 30] . In our case the relevant function is 
— fc^ + C"^ (jf) / rf \ we have a turning point which will de- 
pend on the explicit form of g{'q) and a pole of order 2 
in the limit 77 ^ 0^. Therefore, we have an Airy solu- 
tion around the turning point which goes over to an LG 
solution for conformal time approaching zero. As will be 
shown later, this approximation fits the exact behavior 
of the solution accurately. For the two approximating 
solutions (Airy and LG), the convergence criteria are es- 
tablished in the following way. Assume that we have a 
pole of order n > 2 at a finite endpoint 02 and that g{rf) 
and 9(77) are meromorphic functions of the form 






(02 - vf jr'. 



(12) 



Then the error control function (see Appendix 0) con- 
verges if go = —1/4. For a detailed proof of this fact, see 
Ref. [23 • For the error control criteria to be satisfied in 
our case, we must make the choices 
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- "s 



liv) 



4M 

j^2 
1 



(13) 
(14) 



The final form of Eqn. Q then becomes 



77^ 



'M 



Uk, 



(15) 



where Vg — [z" / z)!^ + 1/4 and the turning point is at 
^"^ — '^sivs)/Vs- Note that the turning point, 775', is a 
given function of k. 

In an exactly analogous fashion, the equation for the 
tensor modes ^ is written in the form 



4(??) - 



Vk, 



(16) 



where i^^ — {a" /a)rf' + 1/4 and the turning point is at 



'■^{vt)/V[ 



''t- 



Unlike the approach of matching solutions through re- 
gions I, II, and III, the idea of a uniform approximation 
is to provide a single approximating solution which con- 
verges uniformly in all three regions with a global, fi- 
nite error bound. The normalization is determined once 
and for all by simply matching to the exact solution as 

k —^ CXD. 

To continue, we follow Olver in defining a new indepen- 
dent variable ^ and a new dependent variable U, given 
by0 



\dr] J ' ' Xdv 



-1/2 



U. 



In terms of the new variables, Eqn. (|ll|l becomes 



d^s, 



[b\ + ^{C)\U. 



(17) 



(18) 



where 



V'(C) = Mv)g"{v)-^g'\v)\ 
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16ff3(r;) 



5(77) 16^2 



-^3/2 ^ 

3^ 



\/9{v)dv, V>V 



i-0 
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y/~g{r])dr], 7; < 77. 



(19) 
(20) 
(21) 



Now imagine neglecting "0(0 ^^ ^ S'^sl approximation; 
then the solution to the differential equation is given im- 
mediately in terms of Airy functions. In the next or- 
der where tp{^) is no longer neglected, the derivation of 
the solution u becomes more involved. Fortunately, in 
Ref. [231 the general solution for Eqn. (|11|) in the uni- 
form approximation to all orders is derived with error 
bounds to be 
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s=0 



(2) 
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i(^0 



giv) 



54/3 

-1-1/4 



C 



Bi(62/3^)f]Ml 



s=0 



Bi'(62/3^)'gS,(e) 
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t2n+l 



(23) 
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s=0 

with coefficients defined by an iterative procedure 

^o(0 = 
BsiO = 



1 w.l.o.g., 



2V±|J0 



±v 



(24) 
,(25) 

(26) 



where the upper signs are to be taken on the right of 
the turning point 77 and the lower signs on the left of 

the turning point. The error terms i2n+i ^'^"^ ^2n+i ^^^ 
discussed in Appendix In the next and following Sec- 
tions, we will calculate the leading and next-to- leading 
order solution for Eqn. (|15|l with explicit error bounds 
and derive the corresponding power spectra and spectral 
indices. 



IV. RESULTS AT LEADING ORDER 

We now turn to the specific form of the approximating 
solutions at leading order. Taking n = in Eqns. 1)22(1 
and H23() we find a solution for Ukiij) containing a part 
valid to the left of the turning point (rj < fj) and a part 
valid to the right of the turning point (rj > fj). The 
unnormalized solutions are 
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y'kljv) = [/s(fc,7y)/.9s(fc,7?)]'/'Bi[/s(fc,,7)], (29) 

^ 2/3 



9s{k,v) 



2 






ns 



V 

-k\ 



j^Z 



^fc.i.< 



1 



(^7)1 < T 



MJU) 



{' 



AV,,c(£) 



'} 



,(30) 

(31) 
(32) 



A i?(/s)Ai(/s=) 

i4i,('?)i < r^tii/y^ -""^^'^-^}' ^^^) 

where the lower index 1 reflects the order of the approxi- 
mation, the functions with index < are taken on the left 



of the turning point, and those with the index > are to 
be taken on the right of the turning point. Mirf), Nirf), 
Va.f3, and A are defined in Appendix ^ and the error 
control function £(77) is given by 
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dv i 9 s 
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(34) 



Inserting the explicit expression for gs(k,r]), Eqn. H31|l. 
and integrating by parts leads to 
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The errors terms ej. j^ ^(77) in Eqns. H27|) and H28|l en- 

capsulate the contributions to u).\l- {if) beyond leading 
order. The general solution for Ukijii) is a linear combina- 
tion of the two fundamental solutions u^ ' {rj) and u], ' (77), 



VIZ., 



- A„W 



Uk{v)^Au'^>{7j)+Bui^>irj) 



(2), 



(36) 



independent of the order of the approximation. In order 
to fix the constants A and B we have to construct a linear 
combination of u^, ' (77) and m^ (77) such that the result 
has the form 7x^(77) — e^''^''/v2fc in the limit k ^ 00. In 
this limit, the domain of interest is region I, far to the 
left of the turning point. In this case, for well-behaved 
'^Si f<{k, 77) is large and negative and we can employ the 
asymptotic forms 
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Making the choices, 
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e'-, B 



e'4. 



(37) 



(38) 



we find 



■uk,i,<{n) 



lim 

-ki) — >c 
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exp 



^i^[f<{kM''^\, (39) 



which is the required adiabatic form of the solution at 
short wavelengths and as 77 ^ —00 |21|. C is a con- 
stant phase factor which is irrelevant when computing 
the power spectrum. 

The 77 ^ 0~ limit defines the region of interest for 
calculating power spectra and the associated spectral in- 
dices. In this region, the l/rj^ pole dominates the behav- 
ior of the solutions and the Airy solution goes over to 



the LG solution. The LG form of the solutfon is more 
tractable than the Airy form, leading to simple expres- 
sions for the spectral indices. We now demonstrate how 
the Airy solution for small 77 approaches the LG solu- 
tion. The linear combination of Eqns. H27I) and (|28l) in 
first order with the appropriate normalization is given by 

«fc4.sW = ]llcf'J\k,rj)g-'/\k,rj) [Ai(/s) - zBi(/j,)] , 

(40) 
with the error bound 

\ek,i,dv)\ < ^{[cxp(AV„„ (£))-!] 
+ [exp(AV;3^„(£)) - 1]} 

derived from Eqns. (EH, I®, IHI), and ^^. (The ex- 
plicit form of the variation of £ will be discussed below.) 
For small rj we are on the right of the turning point; the 
argument of the Airy functions, i.e., />(fc,?y), becomes 
large and the Airy functions can be approximated by 
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which leads to 
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(43) 



For computing the power spectra in the krj ^ limit, 
only the growing part of the solution is relevant: 



Uk,i,>{v) = 1™ -iC 



Ort-HIl^^"-""''' 



(44) 

Once the approximate solutions to Eqns. Q and (Q 
have been found in the manner described above, the rel- 
evant power spectra can easily be computed. The def- 
inition of the scalar power spectrum for C, where it is 
understood that all time-dependent quantities are to be 
evaluated in the limit 77 — > 0~, is 
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Pi,s{k) denoting the power spectrum for the scalar per- 
turbations in the leading order approximation and drop- 
ping a second-order term in efe,i,>. Substituting the LG 
expression for Uk from Eqn. (|44l) . we have 



Pi,s{k)= lim 



fc3 



-T] 



kr,^Q- 47r2 |2;(77)|2 1/5(77) 



exp 



[f>ik,v)] 
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(47) 

with the error bound for the power spectrum given in 
Eqn. igni). 

We now discuss the variation of £. (See Appendix IbI 
for a short discussion on the variation of a function in 
general.) We are interested in the error bound of Uk^i 
over the full domain of interest — cx) < 77 < 0~, which 
implies /3 = — cx) and 77 — > 0^. In the general case. 



V-oo,,(f)-5]|f(a)-f(/3)|, 



(48) 



where the sum is over all individual monotonic subintcr- 
vals {a, (3) C (—00, 77) of £. In the special case of mono- 
tonic £ over the full range of 77 the answer can be given in 
a simplified form: By inserting the definition of gs{k,r]) 
from Eqn. (|31|1 into Eqn. H34() and integrating by parts 
we find for the variation of the error control function: 



t^—OQ^Tj \^ ) 




(49) 



where it is understood that 77 has to be taken in the 
limit 77 — > 0^. With this expression the error bound 
for the power spectrum given in Eqn. (|46|l is completely 
determined. 

The calculation for the tensor power spectrum follows 
along the same lines, yielding 



exp 



4 [~ 



f>{k,v) 



3/2 



Pi,T{k) = lim -—r — — - 

(50) 
with the error being controlled exactly in the same man- 
ner as in Eqns. H46(l and (|49|l with the subscripts S ^ T 
and f{k,ri) indicating that gs{k,rj) in Eqn. 129|) has to 
be replaced by grik, rj). 

Before proceeding to the computation of the spectral 
indices, we make a few observations on the nature of the 
error term above. In previous work |27l | we have shown 
that the error term in the case of constant v introduces a 
fc-independent amplitude correction for the power spec- 
trum which does not affect the spectral index; therefore 
the uniform approximation for the spectral index is ex- 
act in that case. In a later section we explicitly discuss 
this calculation. It turns out to be possible to utilize this 
result in a general way, and we turn to this now. 

Suppose that we split the effective potential into two 
terms, writing 



u'ir,) 



1 

4 



1 



^nO)-T + //^(r7)-^^(0). (51) 



Further we define a corresponding splitting of tfie error 
term, ik,i,> = ^o + £: where eg is the error term for the 
case of constant 1^(77), the constant being given by v{0). 
Following Olver, the total error term and the error term 
eo each satisfies an integral equation 29] ; using these we 
easily derive an integral equation for e with a somewhat 
different inhomogeneity. By construction, this inhomo- 
geneity is reduced in size compared to the inhomogeneity 
which appears in the integral equation for the full error 
term e. Using the fact that the error term for constant 
1/(77) is /^-independent, in this way we demonstrate that 
the full error term can be written as a sum of two terms: 
eg which is fc-independent and e which has the property 
that it vanishes identically for constant iy{ri) and satisfies 
an integral equation with a reduced inhomogeneity. Ex- 
plicit calculation, applying the theorem of Olver in the 
slightly generalized context, gives a relative error bound 
of the form 



following form: 



\^iv)\ < 



1 



Bi(/s) 



E(J^)M(J^)^{f^) exp[AV„,,(£)], (52) 



where E[fg^) and M{f^) are the comparison functions 
introduced in the appendix and used in the discussion 
above, and where <&(/§) = Va_^(£ — Sq); here £ is the full 
error control function and £q is the error control function 
for the constant v^rj) approximant. 

This splitting therefore resums the error contributions 
which arise purely from the constant part of 1^(77) and 
separates them into an explicit fc-independent additive 
contribution to the full error, leaving a term which is 
quantitatively smaller and vanishes in the case of con- 
stant 7/(77). This shows why we expect the uniform ap- 
proximation to the spectral index to be a much better ap- 
proximation than a direct application of the Olver theory 
would indicate. 

Next we discuss the evaluation of the spectral indices. 
The generalized spectral index for scalar perturbations 
can be obtained from the power spectrum via 



ns{k) = l + 



dlnPs 
dink 



(53) 



Differentiation of the power spectrum with respect to k is 
straightforward. It is important to note that the turning 
point fjs is a function of fc, since fc — ~i's{vs)/\ils\ where 
vsivs) is the value of i^siv) ^t the turning point t] = fjs- 
Using this relation, one finds 



nis{k) = A-2k^ lim 



dr]' 



vs V9s{k,v') 



(54) 



Following from the discussion above, the error in the 
spectral index arises only from the fc-dependent part of 
the error in the power spectrum, which vanishes in the 
case of constant vs- Thus the error in the spectral index 
is sensitive only to the time variation of vs- To estimate 
this error, the spectral index as written in Eqn. (|53|1 can 
be expressed via the leading order power spectrum in the 



nsik) = 1 + 



dHPi,S + e^.^s) 



1 



dlnfc 



d\nk 

?^i,s(fc) + efe,i,s> 



d^k,i,s 
dk 



with 



de 



^k,l.,S 



k,l,S 



dk 



(55) 



(56) 



Here we have neglected error terms of order e^. We can 
estimate je^ 1 5I by using the next-order contribution to 
the power spectrum (Section [VJ, leading to the final re- 
sult 

fc2 



^liMn) 



r/dr] 

V27 (4 - fc2772)3/2 



(57) 



i-5!yi 



Vg — k'^rf 



Again, it is understood that 77 has to be taken in the limit 
7; -^ 0^, thus the spectral index defined in Eqn. (|54|) has 
no time-dependence. 

The above analysis can be carried out for tensor per- 
turbations in an identical fashion, including the error es- 
timation, with the replacement vs —^ vt- The spectral 
index for gravitational waves is given by 

di 



rv 
ni.T(fc) = 3-2fc2 lim / 



, , (58) 

'f,T V9T{k,f]') 

Approximate evaluation of the integrals in Eqns. I|54|) 
and (|58|l can be performed to yield more familiar forms 
for the spectral indices as shown in Section IVll 



V. NEXT-TO-LEADING ORDER 

To proceed to the next order, we first write the unnor- 
malized solution for Uk, following directly from Eqns. H22|) 
and idSl: 

^k.Uv) = [Uik,v)/9s{k,v)f^{Ai[U{Kr^)]{Ao[U{k,fj)] 

+ Ai[/5;(fc,77)])+Ai'[/s(fc,ry)]i?o[/s(fc,^)]}, (59) 

^'sgW = [/r(fc'^)/ff^(fc'^)]'^'{Bi[/s(fc,77)](Ao[/s(fc,7y)] 
+ Ai[/g(fc,77)])+Bi'[/g(fc,7;)]Bo[/g(fc,77)]}, (60) 

with 

Ao[U{k,v)] = 1, (61) 

±1 ffi^ i;{v) 



BoiUik.v)] = 



2v/±/s(fc,77)7o V±v 



zdv. 



1 



AiiUik,!^)] = --B'o[f^{k,7j)] 



(62) 
(63) 



2 / i'[U]Bo[Uik,v)]d[U{k,rj)lBA) 
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The error bounds in next-to-leading order are given by 
[derived from Eqns. (|ijl|l and (|ij2|l ]: 

\ei%\ < 2E-\f^)M{f^)Wf^,p, (65) 

\e^klj < 2E{f^)M{f^)Waj^, (66) 



BiiU) 



mv)A^iv) ^ A'liv)], 



with 

(67) 
W^j^ = exp{2AV„j^(|/s|i/2i3o)}v„j^(|/sr/2^i), 

(68) 

and 

±1 f^S dv 
2^/±f;Jo V^ 

(69) 

A'liv) = -li?-(^;) + i[V;'(^;)i?o(^)+V'(^')i?^(«)], 

(70) 

recursively derived from Eqns. (|25|l and H26|l . 

As in leading order, the general solution Ukirj) is a 
linear combination of uj^ {rf) and u^ (??)• Fortunately, 
we will not have to calculate the normalization again; in 
the limit rj — > — oo the Bunch-Davies vacuum is the exact 
solution of the differential equation for Uk{rj), and, in this 
limit, all corrections from the next-to- leading order terms 
are subdominant and of no interest. 

A further simplification follows from the fact that 

(2) 

only the growing solution, u^, , is relevant to determin- 
ing the power spectrum and spectral index and that we 
can restrict ourselves to the solution for Uk in the limit 
kr] —> 0~ . Employing once again the approximation for 
the Bi-function for large, positive argument, Eqn. H42|l . 
and in addition, the approximation for its derivative. 



Bi'(x) = ^xi/4expf -0:2/3 
a/tt V 3 



the normalized Uk{ri) in the relevant regime is 

'2 



(71) 



Uk,3,>iv) **= ^^C"./- 



-^exp]|[/>(fc,,7)]^/2 



X |l - ^B'oih) + Vf>{k,v)Bo{f>) 

i'if>)Boif>)dif>)\ , (72) 



1 

+ 2 

where the error bound is given by 



f-l/4, 



ek,3,>\ < 2/>^'"exp - [f>ik,Tj)r^ Waj, 



3/2 



(73) 



Analyzing Bq in detail shows that the derivative B'q and 
the integral over tpBo are subdominant in the limit krj — > 



0~. Hence the only term leading to a correction of the 
power spectrum is Bq itself. From the general expression 
(03: the power spectrum at next-to- leading order is 



kri^O 



P2,s(fc) "- P^Mk)\^ + 2Vf>ik,v)Bo[f>ik,v)]\, 



with Pi.s(fc) as defined in Eqn. H47|) and 
1 /"-^^ tp{v) 



Bo{f>) = 



4'iv) 



2V7: 



> Jo 



-dv, 



5 , V {Agsg'^ - 5g'i) 



I6w-^ I65I ^ff'gs 



(74) 

(75) 
(76) 



The first term in Bq [after writing out '4^{v) according to 
Eqn. H76() ] can be integrated immediately. The contribu- 
tion from the lower integration limit, which appears di- 
vergent at a first glance, cancels with contributions from 
the other terms in the integral. This can be shown by 
expanding ^p{v) around zero and integrating explicitly. 
The error bound can be calculated in the same way as in 
leading order. We find 



Ps(fc)- lim —^ 



Uk,3,>{'n) 



z{ri) 



[1 + <3,5W], (77) 



with 



2e 



fe,3,> 



-fc,3,S 



Uk,3,>iv)' 



(78) 



The spectral index at this order, n2^Sj cis calculated from 
its definition H53|) . is given by 



"2,s(fc) = "i,s(fc) + lim 



2k 



diVJ>Bo) 



kn^o- |l + 2v^Bo| 



dk 



(79) 

where ni^s{k) is given by Eqn. H54|) . Evaluating the 
derivative of Bq with respect to k leads to the follow- 
ing expression for the spectral index: 



«2,s(fc) == "-l,s(fc) + 



2fc2 



2|l + 2V77Bo| 



(80) 



9s 

/5s V2.9s 
dv 



1 



ffl 



4u2 



-0(w) 



15 

8^ 



Arfgs 

vs 



The error estimate for the spectral index can be obtained 
in a similar way as in leading order, just as for the power 
spectrum. We find 



with 



ns{k) ^n2,s{k) + el^^^s{ri), 
d ek,3,>iv) 



4,3Av) = 2fc 



dkuk^3,>{v)' 



(81) 



(82) 



In the case of the errors in next-to-leading order we have 
to evaluate Waj^^ which is defined in Eqn. H68I) . This 
can be done in principle but the result is rather long and 
complicated and we do not write it out here explicitly. 
In a forthcoming paper |32| we will examine different in- 
flation models numerically and show the corresponding 
results for the next-to-leading order error estimates. 

Proceeding in the same way as for the scalar perturba- 
tions the power spectrum and the spectral index for the 
tensor perturbations, including error bounds, can be cal- 
culated. P2,T(fc) can be obtained from Eqns. H74|) - (|76|l by 
replacing Pi.s{k) with Pi.rik) on the r.h.s. of Eqn. H74() 
and replacing gs{k,r]) and its derivatives in Eqn. H76|l 
by gT{k,r]) and its derivatives. The spectral index for 
the tensor perturbations can easily be derived by replac- 
ing ni^s{k) and gs{k,ri) and its derivatives in Eqn. H8U|) 
on the right hand side by ni^xik) and gT{k,ri) and its 
derivatives. 



VI. LOCAL APPROXIMATIONS 

A. Uniform Approximation 

The solutions obtained so far from the uniform ap- 
proximation are nonlocal; for the sake of simplicity and 
in order to compare with conventional slow- roll results, 
local expressions are desirable even if some accuracy is 
sacrificed thereby. This requires making additional ap- 
proximations regarding the integrals in Eqns. H54I) and 
(|58|l . The analysis is again identical for the scalar and 
tensor cases so we address the scalar case first. The inte- 
grand has a square-root singularity at the turning point, 
i.e., at the lower integral limit. At the upper limit rj 
goes to zero and the integrand vanishes linearly, there- 
fore, assuming vsiv) is well-behaved, we expect the main 
contribution to the integral to arise from the lower limit. 
Combined with the knowledge that i>s is slowly varying 
it is reasonable to expand z/5 around the turning point 
in a Taylor series. To second order in derivatives I'siv) 
reads 



•^Uv) - vl+2vsv'g{iq^ 



fjs) + {i?g+9ps)iv-Vsf, (83) 



where the bar indicates that this quantity has to be eval- 
uated at the turning point. We can now solve the integral 
in Eqn. H54|) exactly and find for the scalar spectral index 



nsik) ~ 4~2i?s\l 



vs 



-Vs 



(-1) 



(84) 



vl 



.{2-7r) + ^il-n) 



^S 



This is a simplification of the leading order result in the 
uniform approximation; at the per cent level of accuracy 
for the spectral index expected at this order, we have 
verified that it is adequate to keep terms up to second 
derivatives of I'siv): ^^ i^^ Eqn. l|SS|) . 



For the tensor spectral index we find analogously 
nrik) ~ 3-2z^t|i-^^t(i-|) (85) 

2 






B. SloAV-Roll and its Variants 

As a prelude to the comparison of slow-roll results with 
those obtained from the local approximation discussed 
above, we give a brief overview of the slow-roll paradigm 
for calculating the power spectrum. Returning to the 
matching problem discussed in Section II (represented 
by Fig.^, a rough estimate of the power spectrum may 
be obtained by extrapolating the high and low-frequency 
solutions to an intermediate regime — ^77 = 1 or fc = 
aH ('horizon-crossing' for the fc-mode of interest) and 
equating them at that point. It is understood that \Ak\ is 
determined from the matching condition at /c = aH , i.e., 
I^fel = l/(2v2fc). With this substitution, the familiar 
leading-order result [33| 



^'»'-lS)(? 



(86) 



k=aH 



is obtained. This expression is useful in providing quick 
estimates for the power spectrum but is clearly not very 
precise: the solutions are being matched in a region where 
they were not meant to be applied. The conventional 
'slow-roir approach to proceeding further is to improve 
the matching by providing a better intermediate solution 
in the region —krj^ 1 or fc ~ aH (region II of Fig.^ and 
then to match the short and long wavelength solutions 
against the intermediate solution. 

To understand the situation more concretely, it is use- 
ful to express z" / z and a" /a in the exact forms 



2a^H^ [ 1 + e + ^<5i + 2e6i + e^ + ^62 ) ,(87) 



where 



2a'H^l-\e], 



- Nt , 6„ 




1 d" 



H^d) dt" 



(89) 



Here we have followed the notation of Stewart and 
Gong |34| . which is especially convenient for comparing 
results expanded order by order in slow-roll parameters 
as given in the next section. (This convention is slightly 
different from that used in our previous study 27] which 
followed the conventions of Ref. |13|.') 
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The derivatives of the leading order slow-roll parame- 
ters e and ^i are of second order in the slow-roll param- 
eters: 



H 



2Sie + 2e^, 
62 + 616 -Si. 



(90) 
(91) 



In principle, the parameters e and S„ are functions of 
time that can be determined by solving the Friedmann 
equations. In the slow-roll approximation, the values of 
e and 6n are assumed to be small, and one aims to solve 
Eqns. (0J and (O in terms of expansions in these param- 
eters. At leading order, with e and Si both -C 1, S2 being 
already next order can be neglected: 62 ~ 0(e^,(5i, e(5i), 
then it follows immediately from Eqns. H90|l and H91|) 
that, to leading order, the derivatives of e and i5i are 
approximately zero, and therefore they can be treated as 
constants. It is crucial to note that at higher order in 
the slow-roll expansion, the assumption that the param- 
eters are approximately constant does not hold, and the 
analysis becomes complicated [23|. 

With these results in hand, the leading order slow-roll 
versions of Eqns. I@J and become 



< + 


e- 


1 
-2 


(.- 


IV 

4J. 


Uk = 


= 0, 


(92) 


v'^ + 


e- 


1 


\ 


-i): 


Vk = 


= 0, 


(93) 



where the leading order vs and vt are given by 

3 3 

vs = - + 2e + 6i, VT = -+€-. 



(94) 



These equations can be solved in terms of Hankel func- 
tions as 

Mv) = ^e^(''^+3)f (_^)iiJ^^(i)(_fc,?), (95) 

vk{v) = ^e^t-^^+s)! (_^)^i/^^(i)(_fc^). (96) 

The Hankel function solutions can now be matched to 
the short and long wavelength solutions as described in 
Ref. [23 (see also Ref. 24] ) . To calculate the power spec- 
trum we need in addition an expression that connects 
the conformal time rj with the slow-roll parameters. This 
can be obtained to any order by repeated integration by 
parts p^: 



aH L 
+2{)e^5i 



£ + 36^ + 2eSi 
2eSl + 2eS2 + ■ 



Ibe-" 



(97) 



It follows from Eqn. H97|) that at leading order, 77 ~ 
— (1 + e)/aH . Using this result and the small argument 
approximation for the Hankel functions, the first-order 



slow- roll corrections to the scalar power spectrum (|86|l . 
as originally computed by Stewart and Lyth J33, read: 



Ps{k)~[l-{2C+l)€ + CSi\ 



H^ 



k=aH 



(98) 

where C ~ —2 + ln2 + 7 ~ —0.73 and 7 is Euler's con- 
stant. Errors inherent to the matching procedure and use 
of the Hankel approximation are discussed in Refs. |23| 
and H. 

As already explained above, straightforward extension 
of the slow-roll app roach to higher orders runs into seri- 
ous difficulties |23l |24| . The primary obstacle is the fact 
that at next-to-leading order, as is clear from Eqn. H90|l . 
the slow-roll parameters can no longer be treated as con- 
stants and the Hankel solution no longer holds. Thus, 
to obtain more accurate results, one has to abandon the 
original slow-roll expansion technique of fixing a Hankel 
solution in the intermediate region and aiming to obtain 
higher-order expressions for the Hankel index in terms of 
an expansion in slow-roll parameters. 

An alternative approach within the slow-roll method- 
ology is the work of Stewart and Gong |^ . In this calcu- 
lational scheme, explicit matching of solutions is avoided 
by using a perturbative approximation, allowing the in- 
clusion of higher-order corrections (The leading order so- 
lution is still the Bessel approximation of Ref. 35].) To 
be more concrete, Stewart and Gong employ the following 
approach: instead of replacing z" / z in Eqn. Q by a con- 
stant divided by 77^ as in the original slow-approximation, 
they choose the ansatz z — t]^^ f{hir]). This ansatz pos- 
sesses, in addition to the term proportional to 1/77^, a 
contribution which depends logarithmically on time. In 
the differential equation for Uk this additional term is 
then treated as an inhomogeneity and the solution for Uk 
can be found using Green's function methods. The final 
expression for Uk is an integral equation. The integral in 
the expression for Uk is then expanded in slow-roll pa- 
rameters. By including a weak time dependence in z" /z 
this approach allows the extension of the slow-roll ap- 
proximation to higher orders without violating Eqn. I|90|) . 
Nevertheless, this approach also lacks an error estimate. 
In a second paper [23, Stewart elaborates on these re- 
sults, though only at leading order, by choosing a gen- 
eral expansion point for the slow-roll parameters (instead 
of fc = aH) and showing that expressing the slow-roll 
parameters in terms of potential derivatives (see, e.g., 
Ref. [33 for a discussion of this approach) can lead to 
incorrect results. An extension of this work to one more 
order has been recently performed |37| . 



C. Slow-Roll Redux 

We are now in a position to discuss our results in 
the context of the slow-roll analysis of the previous sub- 
section. In our case, obtaining spectral indices from 
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Eqns. H84() and H85() as a function of e and (5„ is straight- 
forward: we simply expand vs, vt, and their derivatives 
in terms of slow-roll parameters. A crucial point to keep 
in mind is the value of k in terms of which the results 
are stated. In the slow-roll expansion the evaluation of 
the slow-roll parameters is traditionally given at horizon 
crossing, —ki] ~ 1, a somewhat arbitrary choice due to 
the uncontrolled nature of the approximation. As shown 
in detail in Ref. [23 calculating the slow-roll parameters 
at a convenient time close to horizon crossing leads to 
small finite corrections in power spectra and spectral in- 
dices. In our case, the natural expansion point is the 
turning point, therefore a truly direct comparison with 
slow- roll results, such as those of Ref. PJ]) is rather com- 
plicated. 

We begin by considering the spectral index for the 
scalar perturbations as given in Eqn. (JHIJ. In order to 
write ns in slow-roll parameters up to second order we 
expand vs, I'sVs, f^nd Vgfjg up to third order in slow- 

roll parameters. The contribution Vg fjsl^s is already of 
fourth order and is neglected. Using the expression for 
conformal time in slow- roll parameters given in Eqn. I|97|l . 
the relation between ^^5 and z" / z, 



335 



9 Z 9 1 

^^ = T^ +4' 



(99) 



and the expression for z" / z given in Eqn. (|87|l . we find 
for the different contributions: 

3 o ? 16 2 14 - 1-2 
1- 206 3 2-3 -2 287 2? 



2fi - 2 - - 



(100) 



r^sv's = -46- - 56-^1 +51-5^- ^^iH, ~ ^e^ - ml 

19 - 2 - - - 1 - 
-ye^2- 3^? +^1^2 -3^3 + 0(6*), (101) 

fsv'l; = M^ + hllx - 5? + (52 + — -e^ + ^-e^(5i - \Exh 

+ ^eJ2 + h\ + m\ -K ^^3 + 0{^). (102) 

As in earlier sections, the bar indicates that the slow-roll 
parameters are to be calculated at the turning point, and 
0(t^) represents all slow-roll terms of fourth order and 
higher. Finally, inserting Eqns. H1UU|I - H1U2|) into the local 
expression for the scalar spectral index (|84|) allows us to 
write the spectral index in terms of slow-roll parameters: 



ns{k) ~ l-4e-2^i-8e2 



-mig-. 



17 



2(<5f - <52) ( y - ^ 
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3 I 15 / 3 H3 



iHr 



1942 



1005 



- mil 



20 



59 _^ /286 
3 ^ I5 



5^1(52 ( 4^ - ^ 



(103) 



Thus, in the end, starting from the nonlocal expression 
for the scalar spectral index given by the uniform ap- 
proximation [Eqn. (|54|l ]. we have finally arrived at a local 
expression for ns in terms of slow-roll parameters by em- 
ploying two expansions: First we expanded the integrand 
in the expression for the spectral index in a derivative ex- 
pansion in vs, to solve the integral in Eqn. H54|l . Then 
we further expanded the result in slow-roll parameters. 
However, the two expansions arc not independent; had 
we decided to stop the expansion in derivatives in ^^5 af- 
ter the first term we would not have obtained the second 
order slow-roll contributions from the expansion of the 
second derivative of vs- If one wants results quoted to 
some order in slow-roll parameters, this requires going 
up to a finite order in derivatives of vs', however, written 
in this way, it is also clear that a 'priori it is not obvious 
which expansion is the dominant one - the expansion in 
derivatives of vs or the expansion in slow-roll parameters. 
In the absence of further information regarding vs it- 
self, the question cannot be answered satisfactorily. This 
demonstrates one of the inherent difficulties of deriving 
higher order expressions for the spectral index via Taylor 
expansions without having a well-defined error bound. 

The analogous result in Ref. |3J| agrees at leading or- 
der. The forms of the higher order contributions are 
apparently different due to the difference in evaluation 
points and the different approximations employed. How- 
ever, both results should be treated with some caution: 
(i) Without an error control theory, it is not clear that 
inclusion of higher order terms actually improves the ac- 
curacy of the result (convergence is not guaranteed since 
the Taylor expansion leads only to an asymptotic expan- 
sion) (ii) the evaluation point of the slow-roll parameters 
leads to an uncertainty in the calculation - if one really 
wants results accurate to the per cent level, this uncer- 
tainty is important. In order to obtain results with high 
accuracy, error controlled approximations appear to be 
necessary. 

For completeness we also give the result for the tensor 
spectral index expressed in slow-roll parameters. Since 
Stewart and Gong do not derive an equivalent expression, 
we compare it instead to the quadratic-order slow-roll re- 
sult obtained originally by Stewart and Lyth |3a] . For the 
tensor spectral index derived from a slow-roll expansion 
of the local result Eqn. (|85|l . we find 



nrik) 



-2e-2 



23 



14 



n]eSi, (104) 



whereas Ref. [35l| obtains 

nrik) = -2e - (3 + c)e^ + (1 + c)e6i 



(105) 
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with c — 0.08145. As for the scalar spectral index, the 
leading order contributions agree. 



VII. THE SPECIAL CASE OF CONSTANT v 

In order to demonstrate the accuracy of our approx- 
imation explicitly we now investigate a special class of 
exactly solvable inflation models, where vg and vt are 
constant. This class includes models such as power-law 
inflation or inflation near a potential maximum. We re- 
strict ourselves in this section to scalar perturbations; 
tensor perturbations can be treated in the same way. The 
exact power spectrum for time-independent vg evaluated 
from the general expression (|45|l is 



Pf(k) 



T;^Y\vs){-A i-kr.r'-'sk^ (106) 
2n-^ \a(j)/ 



X 1 + 



1 



1 



%vs 72v1, 



(107) 



where we have used Stirling's formula to replace the T 
function. [For a detailed derivation of Eqn. (|106f) see, 
e.g., Ref. [l3-] The spectral index is easily found from 
Eqn. (jSnj: 



ns = 4- 2vs. 



(108) 



The general expression for the power spectrum in leading 
and next-to-leading order in the uniform approximation 
is given by Eqn. (|47|l and Eqn. (|74|l . respectively. For 
constant vsi the integrals which appear in these expres- 
sions can be solved exactly, leading to the results: 



P'i\k) = 



221/5-2 



-2I.S 21/s-l 

e Vg 



,(2), 



P'^'ik) 



22^3- 



-e-^-^yf^-^ 



1 

&vg 



H 

a(f> 

H 

acf) 



{-kr^f-^^^e, 
(109) 

2 

(-fcr;)i-2-sfc2 
(110) 



Comparison of these results with the exact power spec- 
trum reveals a nice feature of our approximation for this 
special case: Improving our approximation order by or- 
der leads to matching corrections to the F-function in 
powers of inverse v. 

A more rigorous analysis of the errors in the uniform 
approximation along the lines explained in Sections IIVI 
and FVl shows that the leading order solution is bounded 
by the absolute value of the relative error 



£1,2 



< %/2 



1 

Qvg 



A 



724 



(111) 



where A ~ 1.04. The error in the power spectrum given 
in Eqn. (|109|1 falls comfortably within the bound. 

Similarly, we can solve the integrals in the expressions 
for the spectral indices in leading and next-to-leading 
order, Eqns. lfM|) and (|5n|l . exactly. We find that the 
spectral index in leading order is already exact: 



,(1) 



,(2) 



^A-2vg. 



(112) 



This is an important result demonstrating the high ac- 
curacy of the uniform approximation already at leading 
order. There are no higher order corrections, to be ex- 
pected since the corrections to the power spectrum in 
second order, Eqn. 1110|l . are fc- independent. We can 
also evaluate the error bound for the spectral index from 
Eqn. H56|l for constant v. Consistent with obtaining the 
exact value of ng, we find that in the small k-q limit, the 
error vanishes. 

In order to demonstrate quantitatively the high ac- 
curacy of the uniform approximation, we now consider 
the case of power-law inflation, where v is constant. In 
this case the scale factor evolves as a power-law in time, 
a (X tP, hence the Hubble parameter 



H{t) = pt- 



(113) 



The power p and vg are connected via 



J^S = 



p-1 



(114) 



The conformal time is given by the exact expression 



V 



1 



1 



aH{l/p^l) 



(115) 



Slow-roll is known to be inaccurate for small values of 
p, we therefore pick p = 2 as a test case - equivalent 
to 1^5 = 5/2. Expressions for the exact power spectrum 
and the ones obtained with the uniform approximation 
are already given above; we now compare these with the 
conventional slow-roll results. Traditionally results for 
slow-roll inflation are calculated at fc = aH. Following 
this choice, the result of Stewart and Gong 34\ for power- 
law inflation is 



P|«(fc) 



H^ 



47r202 
2 

T2 



P 



-(c+1) 



p^ 



TT 





C 1 

2 4 



(116) 



Figs. |2] and O show a comparison of the exact power 
spectrum versus wave number, with leading and next-to- 
leading approximations from the slow-roll and the uni- 
form approximation. The solid black line in both Fig- 
ures represents the exact solution. In Fig. |21 we com- 
pare the first order approximations to the exact solution. 
The dotted line shows the slow- roll result, which has an 
error of almost 30%, while the dashed line shows the 



13 



power spectrum obtained from the uniform approxima- 
tion with an error of roughly 7%. Fig. O displays the 
next-to-leading order result. The slow-roll result (dotted 
line) has improved but is still inaccurate to 9%, which 
is worse than the result from the uniform approximation 
at leading order. The grey dashed line (almost indistin- 
guishable from the black line) shows the next-to-leading 
order result from the uniform approximation, which has 
the remarkably small error of only 0.4%. This result is 
very encouraging; in a forthcoming paper |32l | we will 
demonstrate the excellent performance of the uniform 
approximation for models which can only be solved nu- 
merically. 



for power law inflation up to second order in slow-roll 
parameters is given by 
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FIG. 2: Comparison of the exact power spectrum for power- 
law inflation with the uniform approximation and the slow- 
roll approximation in leading order. 
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FIG. 3: Comparison of the exact power spectrum for power- 
law inflation with the uniform approximation and the slow- 
roll approximation in next-to-leading order. The uniform ap- 
proximation and the exact solution are indistinguishable in 
this plot. 

In addition we can compare the spectral index for 
power law inflation. As explained earlier, the uniform 
approximation leads to the exact answer for the spectral 
index in the case of constant v. The slow-roll answer 



,5fi 



p p^ 



(117) 



For the case p = 2, the slow-roll approximation clearly 
does not lead to a good answer, as the exact result is 
ns — —1 while the slow-roll answer is ns = —0.5. Of 
course p ~ 2 does not represent a realistic model for 
inflation. At higher values of p the slow-roll answer im- 
proves dramatically, nevertheless, unlike the uniform ap- 
proximation, it is never exact. (The amplitude of the 
power spectrum in slow-roll approximation also improves 
rapidly for larger values of p.) 

The remarkable feature of the uniform approximation 
is that the error in the leading order result for the solu- 
tion itself (and therefore for the power spectrum) is al- 
ways ~ 8%, while the next-to-leading order error is less 
than 1%, independent of the choice of p. This robust- 
ness of the error control is a key feature of the uniform 
approximation. 



VIII. COMMENTS ON BESSEL FUNCTION 
APPROXIMANTS 

Before concluding we would like to make another com- 
parison to the commonly used procedure involving Bessel 
functions as approximating solutions for the fundamental 
wave equation. Recall that this approximation arises by 
assuming that i/(r/) is "locally" constant, and it forms the 
basis for methods discussed previously in the literature 

m 

From the standpoint of the Olver uniform theory, it 
turns out that Bessel function approximants arise natu- 
rally when treating the case of an ODE having a coeffi- 
cient function with a pole of order two (29j . Since such a 
pole arises in the model equation at hand, one might ask 
if this theory is not better suited to the computation of 
the power spectrum, perhaps thereby making a connec- 
tion to the methods already known in the literature. 

In fact, direct computation of the first order term in 
the Olver theory with Bessel approximants shows that 
it is not possible to obtain a reliable power spectrum in 
this way. The failure of the method is essentially due 
to the failure of the matching procedure at early times, 
77-^—00. This failure arises because the Bessel approxi- 
mants, although uniformly controlled in ry, lose too much 
of the more delicate fc-dependence of the solution at large 
77, notably for small k. Following Olver's discussion of the 
relation between the Bessel approximants and the Born 
approximation in quantum mechanics [29j |. another way 
to state the problem is to note that the Born approx- 
imation will always fail for sufficiently low energy. In 
essence, Bessel approximants provide a systematic way 
to calculate the power spectrum at high k, whereas we 
need information at low k. 
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The success of the Airy approximants is a resuh of 
a tradeoff. By concentrating on the turning point at 
the expense of some early and late time information, the 
Airy approximants provide a successful bridge for the 
matching operation which is required. One of the costs 
of this tradeoff is a non- vanishing error in the limit r/ ^ 0, 
whereas the Bessel approximants have a vanishing error 
in this limit. However, the error of the Bessel approx- 
imants is not sufhciently controlled in the other limit, 
r] —f —GO, for small k. Furthermore, we have shown 
above how the bulk of the non-vanishing contribution 
to the error for the Airy approximants gives rise to a 
fc-independent amplitude correction, which cannot effect 
the spectral index; this significantly improves a priori ex- 
pectations for the uniform Airy-approximation. 



IX. CONCLUSION 

In this paper we have presented the uniform approx- 
imation as an excellent technique for obtaining power 
spectra and spectral indices from inflation models, mak- 
ing only a minimum number of dynamical assumptions: 
the answers are stated in terms of elementary functions 
and simple integrals which are easy to evaluate numeri- 
cally. The existence of calculable and robust error bounds 
is a crucial advantage of the method. 

In order to completely utilize this approach, the 
next step is a fast numerical implementation, now in 
progress 3Z|. Combined with a numerical code like 
CMBFAST j58j which translates primordial fluctuations 
into the linear radiation and matter power spectra, and 
with our now fully developed error analysis, such a ca- 
pability promises to be very useful for obtaining both 
forward predictions from specific models as well as back- 
ward constraints from observational data. 
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APPENDIX A: GAUGE-INVARIANT 
PERTURBATIONS 

In order to fix notation and to introduce the main 
equation of interest, we summarize results from cosmo- 
logical perturbation theory which are of direct interest 
for inflationary models. The basic formalism is that of 



Bardeen ^; see Ref. [3 for a detailed discussion based 
on this formalism. 

A homogeneous and isotropic Friedmann-Robertson- 
Walker (FRW) universe is described by the metric 

ds — gabdx°'dx 

= a? [f]) [—drf' + "fij {x)dx''dx^] , 



(Al) 



where 7^ is the metric on homogeneous and isotropic spa- 
tial sections, and the conformal time is 77 = J dt' /a{t'). 
Let Di be the covariant derivative on the spatial sec- 
tions. Latin letters denote spatial indices. Considering 
only so-called scalar and tensor perturbations, we write 
the perturbed metric in the form 



ds^ = a{ri)^ [ds's + ds'^] , 
with the scalar and tensor perturbations 



(A2) 



(1 + 2A)dr]^ + 2DiBdr]dx' 


(A3) 


(1 + 2C)7y + 2 (^D,D, - ^7., A<^') E 


dx'dx^ 


jdx'dx^ 


(A4) 



The tensor Eij is transverse, symmetric, and traceless. 
Gauge transformations of this perturbed metric are gen- 
erated by vector fields of the form ^^ = (T, D^L). 

The two standard gauge invariant combinations which 
describe the scalar perturbations of the metric are 



$A = A+{B' + hB)-{E" + hE'), 



$f 



= C--A'''^E + h{B-E'). 



(A5) 
(A6) 



Primes denote differentiation with respect to conformal 
time and h = a' /a. Note that in the longitudinal gauge 
(L — —E, T — B — E'), the perturbed metric becomes 

ds2 ^ q2(^-) ^_(^ _^ 2$4)df?^ + (1 + 2<^chi3dx^dx^] . 

(A7) 
This gauge yields a direct physical interpretation for <i>^ 
and ^c- The tensor perturbations described by Eij are 
themselves gauge-invariant. 

In order to complete the setup for application of the 
Einstein equations, we must introduce matter perturba- 
tions. The energy momentum tensor for an FRW cos- 
mology is that of a homogeneous and isotropic perfect 
fluid 



Tab = Q{'n)UaUb + p{ri){UaUb + gab), 



(A8) 



where the energy density g and the pressure p are func- 
tions only of time and Ua{rj) is the fluid four- velocity, 
which is comoving with the gravitational background, 
Ua ~ —a{ri){dri)a. The most general form for the per- 
turbation of the energy momentum tensor is 

STab = 5gUaUb + 2aqi^aUb)+Sp{UaUb+gab)+Pa^'nab, (A9) 
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where qa is the perturbed velocity field of the fluid, and 
the tensor Hab represents the anisotropies of the pertur- 
bations. The perturbations must satisfy the constraints 



g^'^'qaUb = 0, g'^^Hab - 0, n"^^ = 0. 



(AlO) 



Consider now the case of perturbations with vanishing 
anisotropic stress, Ilab = 0. One consequence of the Ein- 
stein equations in this case is that $ = $yi = ^c- Fur- 
ther assume that the perturbations are adiabatic; then 
the Einstein equations can be combined to give [SJ 

$"+3/i(l+c^)$'-c^A<''$+(2/i'+(l+3c2)(/i2_/C))$ = 0, 

(All) 
where Cs{r])'^ = p' / g' has the interpretation of the adi- 
abatic sound speed squared, and /C is the curvature of 
spatial sections. It turns out that if one defines w by 
$ = 4:TrG{g + pY^'^w, then 



w" - c?A<^'u; - 



-w = 0, 



(A12) 



the only ingredient that we need from this choice is the 
stiff nature of a relativistic field, c^ = 1. We also choose 
flat spatial sections as appropriate for generic inflationary 
models, /C = 0. Since u has a canonical Hamiltonian evo- 
lution, it is a natural choice for canonical quantization; 
the operator is expanded in terms of mode functions as in 
Eqn. (j2Jl, and the mode functions u^ are found to satisfy 
Eqn. Q. 

The tensor perturbations are expanded in terms of ten- 
sorial modes on the spatial sections. 



Etjin,^) 



1 



dpi(k) Wfc(?7)J^ij(x), 



(A18) 



where the mode functions 3^^ (x) have the properties 

(A<^' + e)y,, = 0, 7^^3^y = 0, D^y.j = 0. (A19) 

Again it is found that the mode functions Vk satisfy a 
simple evolution equation, which is given by Eqn. jSj). 



where 



iV^h 



- [h^ -h' + K] 



-1/2 



(A13) 



This is a very convenient evolution equation for the met- 
ric perturbations; it is a complete description of the dy- 
namics of scalar perturbations, under the assumptions 
made on the stress tensor. 

A second form of the evolution equation can be ob- 
tained by introducing a scalar u, which is deflned by the 
relation 



(A14) 



c«(^) l^ + TrlcsW" 



Because of the relation between this definition and the 
differential operator which appears in the w equation, 
it is not difficult to show that u satisfies the evolution 
equation 



z 



where 



"zs (""-'■' + '=> 



1/2 



(A15) 



(A16) 



Note also that these definitions relate the quantity u/ z 
to the metric perturbation in a simple way. 



A'^'>$ = -(47rG) 



1/2. 



h' + 1C (u\' 



hc^^ 



(A17) 



Recalling that $ has the interpretation of the Newtonian 
gravitational potential, this Poisson equation indicates 
how u/ z acts as a source for the potential. 

Finally we specialize to the situation where the back- 
ground is dominated by the dynamics of a scalar field; 



APPENDIX B: DEFINITION OF THE ERRORS 

A key advantage of the uniform approximations pre- 
sented in Ref. [S^ is the uniform control over the remain- 
der terms for the approximations. This uniform control is 
obtained by carefully separating the dominant influences 
in the coefficient functions of the ODE. Because of this 
uniform control, this approach is superior to the earlier 
results of Langer (for a good description see Ref. p?9J). 
Additionally, Olvcr constructs higher order approxima- 
tions not present in the original work. In this appendix 
we review in detail the general error formulae given by 
Olver. The errors in Eqns. I|22|l and (|23|) are bounded by 



i4Vi,g(^e)i |g4»Vi,g(^o/gei 

M(62/3e) ' 62/37V(&2/3^) 



(Bl) 



< 2i.-(6n)exp{?M!^o)\ ^^^^d^l^^^^") 

i41i.g(^e)i |g4lVi,s(^g)/gei 



52n+l 



M(62/3e) ' ^2/3^(^2/3^) 

< 2i.(6n) exp{ ^^^"-^(f^^^°n ^^M\'''B.) 



(B2) 



52n+l 



where M{x) and N{x) are modulus functions, and E{x) 
is a weight function defined as 



M{x) 
M{x) 

N{x) 

N{x) 



= v/2Ai(a;)Bi(a;) for x < c, 
= Y Ai^(a;) + Bi^(x) for x > c, 



(B3) 



.,2 



'2/ 



Ai'^(x)Bi'(a;) + Bi'"(a;)Ai^(x) 



1/2 



Ai(x)Bi(a;) 

1/2 



{Ai''(a;)+Bi''(a:)} 



for X < c, 



for x > c, 



(B4) 
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E{x) 



'Bi(x) 



for c < X < cxD, 



V Ai(x) 
E{x) = 1 for — oo < a; < c, 



(B5) 



and c ~ —0.36605. Some explicit numerical values of 
these functions are given in Ref . {29l] . The auxiliary quan- 
tity A is defined by 



A= sup \7t\x\^/'^M^{x)] 
(-00,00) '- J 



(B6) 



A numerical estimate for A is A ~ 1.04 29]. Finally, in 
Eqns. IjBll) and (|B2|I . we introduced the total variation 
of a function over the interval {a, /3), Va.fiif)- The total 
variation of a function f{x) over an interval (a, (3) is the 
supremum 

ra-l 

Va,/3(/) = sup ^\f{x,+i)- f{Xs)\, 

{a<xo<...<x„<...<P} ^^Q 

(B7) 



for unbounded n and all possible subdivisions of the in- 
terval, a < xq < . . . < x„ < f3. In case of a compact 
interval [a, /?] one possible subdivision is given by n = 1, 
xq = a, and xi — b. Hence 



Va,M/)>|/(/3)-/(a)|. 



(B8) 



Equality holds when f{x) is monotonic over [a, /3]. When 
/(x) is continuously differentiable in [a,/3] we have 



Vc,M/) = / \nx)\dx. 



(B9) 
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